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Abstract 

The resonant transition effects MSW and NSFP for three flavour Majorana 
neutrinos in a supernova are considered, where the transition magnetic moments 
are hkely to play a relevant role in neutrino physics. In this scenario, the deformed 
thermal neutrino distributions are obtained for different choices of the electron-tau 
mixing angle. Detailed predictions for the future large neutrino detectors are also 
given in terms of the ratio between the spectra of recoil electrons for deformed and 
undeformed spectra. 
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1 Introduction 



Though the Standard Model (SM) for fundamental interactions has received a striking 
confirmation from experiments, the neutrino sector still remains a puzzling subject for 
theoretical speculations and for conceiving new experiments. The determination of the 
possible neutrino masses and the consequent mixing angles is the main goal of this huge 
scientific effort. 

Several accelerator experiments have analyzed a quite wide range of the space of neu- 
trino parameters by putting a certain number of upper bounds on these physical quantities 
However, it seems that the most interesting and promising observations should come 
from the cosmological and astrophysical side, which enables to detect extremely weak 
phenomena as neutrino oscillations. 

In this concern, the new experimental results on solar neutrinos seem to confirm the 
previously observed deficit of their flux 0. On this ground, the so-called solar neutrino 
problem can be reasonably considered as a well-established experimental observation. 
This, of course, has a number of implications on the previous question, concerning the 
fundamental parameters of neutrino physics (masses, mixing angles, e.m. properties). 
In particular, the solution of solar neutrino problem needs massive neutrinos and not 
vanishing mixing angles. 

The mechanisms proposed to explain the deficit in the observed number of solar neu- 
trinos with respect to the theoretical predictions of the Standard Solar Model (SSM) 0] 
are basically three. As far as the pure oscillation in vacuum is concerned, it is strongly 
disfavoured by the experimental evidence for a quite strong energy dependence of the 
flux depletion mechanism 0]. Certainly, the most promising solution of the solar neu- 
trino problem is provided by the resonant oscillation mechanism, the Mikheyev-Smirnov- 
Wolfenstein effect (MSW) which has the correct energy dependence. Nevertheless, one 
can imagine the superposition of a Neutrino Spin-Flavour Precession (NSFP) mechanism 
too 0. This effect is relevant when a strong transverse magnetic field is present and if 
the magnetic dipole moments of neutrinos are large enough. 

The nature of neutrinos as Dirac or Majorana particles is also still unknown. In 
scenarios proposed by GUT theories as 5*0(10) neutrinos are described by Majorana 
spinors. In this framework, since new physics is present beyond SM, the e.m. properties 
can also get a strong enhancement. This is in particular the case of some extended unified 
theories 0. 

In this paper we extend the analysis recently performed by some of the present authors 
0. We will consider three-flavour Majorana neutrinos for which, due to CPT invariance, 
only flavour- changing magnetic dipole moments are, in general, non vanishing [||. The 
resonant neutrino oscillations are considered in a supernova where the conditions of large 
density and strong magnetic field are likely to be obtained. 

The paper is organized as follows. In section 2 we consider the effective Hamiltonian 
describing the interactions of three flavour Majorana neutrinos in matter. In particular we 
study the MSW and NSFP effects as a function of the mass difference and mixing angles. 
In section 3, the main supernova characteristics in terms of radial mass density profile 
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and magnetic field are reviewed, and the adiabaticity conditions are studied. The out- 
coming thermal neutrino spectra are then obtained as functions of survival and transition 
probabilities. In section 4 the predictions about deformed thermal neutrino spectra are 
obtained, and the ratio between the spectra of recoil electrons for the deformed neutrino 
distributions and for the unmodified ones is computed. Finally, in section 5 we give our 
conclusions. 



2 Resonant transitions of Majorana neutrinos 

For three fiavour Majorana neutrinos the independent degrees of freedom, which do not 
take mass at the intermediate scale (10^^ GeV) and are coupled with weak interaction, 
can be denoted by , , ^tL and for antineutrinos by , Vfj^n , VrR . For simplicity, 
hereafter we will omit the indication of chirality being clear that it is left-handed for 
neutrinos and right-handed for antineutrinos. 

In presence of a transverse magnetic dipole moment a neutrino can fiip its spin and 
thus change its chirality. However, CPT invariance forbids these transitions for Majorana 
neutrinos, unless they change fiavour at the same time, namely u^l T^fSR, with a,(3 = 
e, fj,, T and a ^ (3. 

This kind of transitions are called Neutrino Spin-Flavour Precessions (NSFP) [^, and 
their probability can be large if the amplitudes of the transverse magnetic field and of the 
off-diagonal dipole magnetic moment are large enough. 

As neutrino oscillations in matter (MSW), even the NSFP can receive a resonant 
enhancement from the presence of a dense medium |]TU|. This is for example the case of 



stellar matter or the extreme condition of a supernova. 

Let us denote with z/ the vector in the fiavour space z/ = (z/g , , i^r ), and analogously 



V = (z/g , z/^ , z/,- ). The 3x3 unitary mixing matrix can be written as |jTT 
U = 



( CipCcj CipSf^ Sip 



We will assume for simplicity 6 = 0, and thus U represents the mixing matrix for antineu- 
trinos as well. 

The evolution equation for neutrino wave functions travelling along the radial coordi- 
nate r ~ ct is 





Ho B±M 
-B^M Ho 




(2) 



The symmetric matrix Ho is the 3x3 hermitian effective Hamiltonian ruling the resonant 
flavour transition in the flavour basis; it is given by 





Ho = -^U I ml \U^ + 




ml 



( Ni{r) 

iV2(r) I , (3) 
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where iVi(r) = y^GpiNe - iVn/2), A^2(r) = -GpNn/V^, and Nn being the electron 
and neutron number density, respectively. Substituting (0) in (H) we get 



[Ho) 



+ N,{t) - KSl - /\C,^Cl , 



-C2ujS2ipC^ 



(Ho] 
(Ho] 
(Ho] 



S + N2[r) 
^ + N2{r)-KClC, 



A 



S2wS2il)S^ + C2c<; ( 5*^5', 



^2 r<2 



^2 



A 



S2u)C2ipStp + -C2u)S2%l) (l + "S", 



2 



A 



(4) 



with S = (m? + ml)/AE, A 



mf)/4E and A = (mf + - 2m|)/4E. The 



corresponding Hamiltonian for antineutrinos, denoted with Hq^ is obtained from (^) by 
replacing Ni[r), N2{r) — > — iVi(r), — N2{r). In what follows the diagonal contribution 
in (^) proportional to S can be neglected since it only gives a common phase factor. 
The quantity M is the matrix of magnetic dipole moments 



/ie/x 
M = I -/ie;. 

The resonant conditions for transformations z/p 




(5) 



and for Up 



> , Uj. can 

be obtained by requiring the approaching of two different eigenvalues of H in Eq.(^. 
For small mixing angles this essentially coincides with the condition of having coincident 
diagonal elements, namely, neglecting second order terms in 5*^,5'^ and S'^ 



iVi(r.)-iV2(r,) =2A 
iVi(r,)-iV2(r,) = A-A 

N^irs) + N2ir,) = 2A 
Ni{r,) + N2irs) = A - A 



(6) 
(7) 
(8) 
(9) 



z/^ , Vr are ob- 



In the same way, the resonant conditions for V^. ^ z7^ , and for z/g 
tained from the charge conjugate transitions by changing sign in the corresponding 

right hand sides. Note that denotes the value of coordinate r for which each resonant 
condition is verified. 

Once that neutrino parameters, namely mixing angles if), and tu, and A and A are 
fixed, one can study the occurrence of resonant conditions and charged conjugate 

ones, by using the particular density profile for electrons and neutrons, contained in A^i (r) 
and N2{r). 

In Ref.[0 the above description for neutrino dynamics in presence of matter has 
been used in order to explain the solar neutrino problem. The authors assume there for 
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simplicity that flavour mixing occurs in the e-/i sector only. This corresponds to take in 
Eq.(|l]) ip^ip « ijj and 5 = 0. Thus, since and Vr are completely decoupled, one can 
consider in only the 2x2 matrix involving e and [i labels. 

To be predictive about neutrino fluxes, one would need the precise expression for the 
solar magnetic fleld. Unfortunately, since this form is quite unknown till now, the authors 
consider several magnetic fleld conflgurations treated in literature. 

As a result of their analysis, it is shown that for a small value of the mixing angle 
Qf.y^ = uj, Ve ^ T^iJ. transitions are sufficient to account for the solar neutrino problem. In 
this case the predictions strongly depend on the magnetic fleld conflgurations even if a best 
flt is achieved for the so-called LIN2 |]12| solar magnetic fleld parametrization. Typical 



values of neutrino parameters able to reproduce the data are Amg^ ~ (10 ^ ^ 10 eV"^ 
and sin (26^,,) < 0.2 ^ 0.3 [111 



3 Thermal neutrino spectra in a supernova 

The neutrino flux dynamics contained in Eq.s(0)-(^) can be also applied to the case of a 
supernova. In a previous paper some of the present authors already considered this 
situation, but in the case of negligible magnetic moments. In the following we extend 
that analysis in order to consider the flavour changing magnetic moment terms too. 

In a supernova the mass density ranges from ~ 10~^ g/cm^ in the external envelope 
up to ~ 10^^ g/cm^ in the dense core. We will assume, hereafter, the following density 
proflle |]TB| 

P-Pof^f , (10) 



r 

with po — 3.5 X 10^^ g/cm^, and Rq ~ 1.02 x 10 ''cm. The quantity Rq denotes the radius 
of the so-called neutrinosphere, which represents the bounding surface of the region in 
which neutrinos of a given flavour are in thermal equilibrium. Note that, the electron 
fraction number Yg outside the neutrinosphere can be assumed almost constant and flxed 
at the value Y; = 0.42. 

As far as the transverse magnetic fleld is concerned, one can safely assume the simple 
expression [|T^ 

i?x(r)^5of— f , (11) 



r 

where r denotes the radial coordinate, and the constant BqR^ ~ 10^'^Gauss cm^. 

Here, we will assume that inside the neutrinosphere the resonance conditions are not 
satisfled, as suggested by predictions on neutrino mass spectrum of a wide range of unifled 
gauge models. In flrst approximation neutrinos are therefore emitted from this surface 
with a Fermi-Dirac distribution 



0.5546 ^2 



1 + exp ( — 



(12) 
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with the different ffavours equally populated. The index a in Eq.(|l2]) denotes the par- 
ticular neutrino species. Since the production and scattering cross-sections for electron- 
neutrinos are larger than for the other flavours, z/g are produced a bit more copiously 
with respect to the other ones. Thus, their neutrinosphere is larger than that for i/^ , 
Vr . This implies that the temperature Tq, in (0) for the z/g -sphere is lower than the 
one of Ufj^ , . Furthermore, since z/^ and are produced and scatter on the surround- 
ing matter only through neutral currents, they have identical spectra. Obviously, since 
i^e 5 ^'/x 5 and Ve , , Vr are produced in pairs, the magnitude of neutrino and an- 
tineutrino distributions are equal for each flavour. For the temperature of z/g and z/^ , 

neutrinosphere we adopt the typical values Tg ~ 3 MeV and = ~ 6 MeV . As 
far as the z/g distribution is concerned, it is characterized by Tg ~ 4 MeV, whereas for 
z7^ and we have = Tf ~ 6 MeV. 

Assuming that the solar neutrino problem is solved in terms of z/g NSFP, as 

explained in Ref. |T2[, from Eq.s(|^)-(|^) we see that in a supernova (for Yg = 0.42 outside 
the neutrinosphere) four resonance conditions can be fulfilled. With decreasing density, 
and for r neutrinos more massive than fi ones, we first encounter the region for e-r 
resonance transitions, and then that for e-/i ones. In order we have: z7g z/^ , z/g ^ , 
T^e ^ i^fi and z/g ^ . More in detail, in the following diagram, the MSW and NSFP 
transitions occurring in a supernova are reported. 



> Ur 

- > z/„ 



(13) 



ly,, 



> Z/p 



> Up 



> z/g 



- > l^tl 

- > Z7g 

- > Z/g 

> Ut 

- > Z7p 



> z/„ 



> Z/g 



- > Z/g 
> l^e 

- > l^u 



(14) 



(15) 



z/g 



> T^e 

> Vfi 

- > Z/g 

- > Vr 
-> 



(16) 



To deduce the total transition probabilities, it is important to establish if the different 
resonance regions overlap, namely if the resonance widths are larger than their separation 
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in r. In many GUT models it is natural to expect m//^ << rrii/^ . In this paper we make 
this assumption, so that the resonance regions involving e-r flavours and those involving 
e-fi flavours are well separated between them. For example, for 10 MeV neutrinos with 
mu^ — 10~^ eV and mu^ ~ 10 eV we have the e-r resonances around the region of 
density ~ 10^ g/cm? (deep in the supernova); on the contrary, the resonances e-/i is in the 
external envelope of supernova {k. 1 g/cm?). Further, the MSW-NSFP resonance regions 
non-overlapping condition is given by H 



Lp |tan26'Q^| + 



{Amlp/2E) 



Lp<r2- ri 



(17) 



where ri is the radial position of the NSFP resonance, while r2 is that of the MSW one. 



and Lp = — [{1/ p)dp/dry^ is assumed to vary slowly between ri and r2- In Eq.(|T 
and P label the flavours involved in the resonance (a,/3 = e,fi,T), 6ap denote the mixing 
angles, and Am^^ is the squared mass differences of the relevant mass eigenstates. For 
supernova neutrinos (-E ~ ^ 50 MeV), assuming (|I0|) and (|ll]), it is easy to see that 
for all models considered below, the non-overlapping condition (^) is well satisfled for 



a 



both Ve , Ve 



p-r and Z/e , Ue 



Up transitions. Thus, we can conclude that the four 



resonances occurring in a supernova are quite well separated one each other, and as we 
shall see, this implies a considerable simpliflcation, since each resonance can be considered 
independently from the others. 

For the MSW resonance we have a simple semi-analytical formula for the survival 
probability Jll], |16 



a/3 ' 



1 ^1 

2 + 12 



exp 



TT 



af3 



COS 26 COS 29\ 



(18) 



where the adiabaticity parameter 7q,/3 is given by 



lap 



Am^^Lp sin^ 29, 



2Ecos29a^ 



res 



(19) 



evaluated at the resonance point. In Eq.(|T8|) is the effective mixing angle in matter 
P, H, |l6l , and is given by 



tan2C^ 



2Am^^ sin 29^^ 



2V2EGFNe - Amis cos 29, 



Tap 



(20) 



For NSFP transitions [16 

-Pa/3 ij^a T^a) 



1 

2 + 



2 -exp 



TT 

■^lapEalS 



COS 29af3 COS 26*! 



al3 



lap 



SEfilpBlLp 



Am2^ 



res 



(21) 
(22) 
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In this case the effective mixing angle in matter is 



tan 261 



(5/3 



2V2EGf {Ne - Nn) - ^mlp cos 29, 



(23) 



Note that for both the adiabatic parameters 70/3 and 7^/3 of Eg. ([T9|) and (^2]), the adiabatic 
condition ^ ^] reads 

7a/3,7a/3>>l • (24) 

The non adiabatic correction factor of Eq.s([T8|) and (plj), neglecting non adiabatic 
effects induced by the magnetic field with respect to those due to density 0, is given by 
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F{e^p) ~ (1 - tan2(^„^ 
+ 



log(l 



tan 



(25) 



In terms of survival or transition probabilities, referring to diagrams (pr^)-([T6|), it is pos- 
sible to write the expressions for outcoming neutrino distributions as follows 



Ve nr 



nu, +n,^ = [1 - P (z/e ^ z/e )] < +[1-P(17, 

+ [P (z/e Z/e ) + P {Ve ^ )] , 
n^, = P(Z7,^Z7,X +[1-P(Z7,^Z7,)]< 



+ [l-P(z/e--Z/e)]< +P(z7,^Z/e)<^ ,(26) 



v/) - P f Z/p 



^e)]< 



(27) 
(28) 



where n?. 



n° . Note that to obtain the above expression we have only used 



the unitarity and the observation that P (z/g z7e ) is vanishing (up to the first order in 
the mixing angle). 

According to diagrams ([T3|)-(|16D the survival probabilities P {ve ) and P (z/g 
can be written as products of the single survival probabilities at the resonances, namely 



P{ye 
P{Ve 



Per [Ve 
PpT {Vp 



Vp 



Vp 



) PeM {pe 



• Vp 



(29) 
(30) 



while the transition probability P ivf, Vf. ) takes the expression 



P{VP 



Vp 



+ Pep. {J^e 



Z/e ) [1 - Pp^ (z/g Up )] [1 - Ppp {I'p )] 
Z/g ) [1 - Pg^ (Z7g Z7g )] [1 - Pg^ (z/g Z/g )] 



(31) 



Note that all the expressions (p9|-(|3lD are obtained under the assumption that the single 
resonances of diagrams (|T3p-(|TB|) are well separated. 
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4 Numerical results 



According to the above results (p6D -(p8D, the deformed thermal neutrino spectra are 
obtained once that neutrino parameters are fixed. 

For the electron-muon sector, the relevant parameters can be fixed according to the 
explanation assumed for the solar neutrino problem. In this paper, since we are interested 
in a possible scenario in which the neutrino electromagnetic properties play the essential 
role we will choose the NSFP explanation fl^ . In this case, the deficit in the solar neutrino 
fiux is mainly due to the conversion z/g , being assumed Majorana neutrinos, which 

is the natural choice occurring in GUT theories where a see-saw mechanism is at work. 
The alternative scenario of a pure MSW explanation was treated by some of the present 
authors in a previous paper 0. 

In the NSFP framework, the values for neutrino parameters able to reproduce the data 
are Am^^ ~ 10"® eV^ and sin (29^^) ~ 0.2 and /i^^ ~ lO^^Vs El- 



Concerning the parameters for the electron-tau sector, they are less constrained. How- 
ever, we can fix Am^^ ~ ml ~ 25 eV"^ in order to be able to identify Ut has the hot 



component of dark matter |jT^, whereas, for the transition magnetic moment we can 



assume /igr to be of the same order of /ig^, since, typically, the enhancement to the elec- 
tromagnetic properties is due to physics beyond the electroweak interactions, which hardly 
distinguishes between r-leptons and /i-leptons. Hence, the only remaining parameters is 
der, for which we will choose three indicative values, namely, 10~^, 10~^ and 10~®. Values 
in this range are for example predicted by SUSY GUT theories where one could 
expect in principle an enhancement of the neutrino electromagnetic properties. 

In Figures 1-3, we show the deformed neutrino distributions for u^, the sum of z/^ and 
z/^ , and Ve distributions, respectively, versus their initial distributions. In all figures, the 
solid line represents the initial distribution and the predictions for the outcoming neutrino 
distributions are obtained for the above three values of Oer- Furthermore, in Figure 4, we 
plot the adiabaticity parameters 'ye^l, Jet and 'Je^l, Ist of eqs. ( p^D and (^21) . In particular 
the dashed line represents 76,- for the two cases Oer = 10"^ (upper line) and O^t- = 10~^ 
(lower line). 

As one can see from Figure 4, for Oer = 10~^ the MSW z/g ^ Ur transition is the only 
one to be adiabatic, and thus the z/g —>■ Vt conversion is the only one which provides 
an efficient reshuffling of initial distributions. It corresponds (dashed line in Figure 1) to 
a depletion of distribution in favour of Vj. (dashed line of Figure 2). Note, however, 
that the adiabaticity of MSW transitions ([l9|) decreases with neutrino energy, whereas it 
increase for NSFP (^). 

The MSW z/g ^ conversion becomes less and less efficient as O^t decreases, as can 
be seen by the other two lines in Figures 1 and 2. For the other values of Oer-, in fact, 
7er << 1. However, as one can see from the dashed-dotted line of Figure 1, corresponding 
to the extremely small value O^t = 10~^, a conversion of Vf. in other kind of neutrinos still 
remains. In fact, in this case the only conversions remaining are the NSFP as one can 
see by observing the antineutrinos spectra of Figure 3. Note that since we have fixed 
A^e^i ~ A^er ~ 10~^^ /i^, all the deformed antineutrino distributions are almost coincident. 



9 



because, in this case, they are almost independent of the on neutrino mixing angles (see 

m 



Eq.s.dlD-dl)). 

The distortions in the neutrino and/or antineutrino energy spectra can be observed 
in the future large neutrino detectors, now under constructions, like SuperKamiokande 
and SNO. To quantify our predictions on directly observable quantities, let us consider 
for example the v-e~ elastic scattering for neutrino detection. A relevant quantity, which 
measures the deformation on neutrino distributions due to resonant transitions is the 
ratio between the recoil electron spectrum for deformed neutrino distributions and the 
undeformed one. It is defined as 



(T, E*^) 



R[T,E, 



th 



(32) 



where 



S (t, ) 









) - H 







0^(-e e- - 



da 



z/e e ) 



d(j 

+ (nv, (E) + riT,^ (E)) — (z7^ (z7, )e- 



dE' 



(33) 



El^ being the neutrino energy threshold, while Emax is the endpoint in the neutrino 
spectra, and (t,EI^^ is defined as (|33D , but for the initial distributions. Note that 
R is independent of the factor k which contains the information about the percentage of 
supernova neutrino reaching the detector and other reduction factors due to the detector 
set up. For the explicit expressions of the differential cross sections occurring in (|33|) see 
for example Ref.[|I^. 

In Figure 5, R is shown for the three values of ^er considered, by assuming for simplicity 
a vanishing neutrino energy threshold and Emax = 50 MeV. 



5 Conclusions 

In this paper we have studied the deformation of thermal neutrino spectra emitted from a 
supernova due to resonant transition effects. In particular we have considered three flavour 
light Majorana neutrinos, as suggested by GUT theories like SO(IO) for the neutrino 
degrees of freedom not taking mass at the intermediate scale (~ lO^^GeV), according to 
the see-saw mechanism 

In the present analysis, which represents an extension of a previous study in which 
MSW transitions were considered only, we have assumed enhanced electromagnetic neu- 
trino properties, as suggested by several unified models and by the NSFP solution to the 
SNP. Thus, the supernova outcoming neutrinos experience spin flavour precessions |Tn] in 
addition to the usual MSW phenomena 0. 
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We have calculated the deformed neutrino spectra once neutrino parameters, like 
masses, mixing angles and transitional magnetic moments are fixed (CPT invariance for- 
bids diagonal magnetic moments for Majorana neutrinos). In this scenario, since we are 
interested in a situation in which the enhanced electromagnetic moments play an essential 



role in neutrino physics, we assume the results obtained in Ref.|]T2[, where the deficit in 
the solar neutrino fiux is explained in terms of the NSFP conversion z/g . Thus, 

from the experimental data one can obtain all the parameters of e-yU sector. As far as the 
e-r sector is concerned, the magnetic transitional moment /Xer can be reasonably fixed 
to be of the order of /Xg^, since the enhancement should be due to physics beyond the 
electroweak interaction and thus it is reasonable to assume that it should contribute in 
a comparable way to all leptonic families. Furthermore, requiring that v^- is the hot 
component of DM [|T^, its mass can be fixed equal to 5 eV . 

For three choices of -v^ mixing angle, we have then explicitly obtained the deformed 
, ^fj. + and Vf> energy spectra, as reported in Figures 1-3. Finally, our predictions 
for a directly observable quantity, such as the ratio R defined in Eq. (|32|) have been given 
in Figure 5. 
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Fig. 1 




Figure 1: The energy spectra for versus energy are reported. The sohd hne represents 
the initial distribution, whereas the dashed hne corresponds to the distorted spectra for 
9eT = 10~^. The dotted hne and dashed-dotted correspond to ^er = 10~^ and 10~^, 
respectively. 



13 



Fig. 2 
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Figure 2: The energy spectra for n^y^ + Hi,^ , with the same notation of Figure 1, 



14 



Fig. 3 




Figure 3: The energy spectra for n-p^ , with the same notation of Figure 1 
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Fig. 4 




Figure 4: Energy dependence of the adiabaticity parameters 7q,/3 and 7^/3 of eqs. (|T9D and 
(P^). The sohd hne represents and the dashed hnes 7er: the upper one corresponds 
to Oer = 10~^ and the lower to O^t = 10^^. The dotted hne is '^ye^l and the dashed-dotted 
represents 7er. 



16 



Fig. 5 
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Figure 5: The ratio R (see Eq.(P^) is plotted for 6(,t = 10"^ (dashed hne), 10~^ (dotted 
hne) and 10~® (dashed-dotted hne). For simphcity a vanishing neutrino energy threshold 
has been used. 
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